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DIVISIBILITY TESTS AND RECURRING DECIMALS IN 
EUCLIDEAN DOMAINS 

APOORVA KHARE 



Abstract. In this article, we try to explain and unify standard divisibility 
tests found in various books. We then look at recurring decimals, and list 
a few of their properties. We show how to compute the number of digits in 
the recurring part of any fraction. Most of these results are accompanied 
by a proof (along with the assumptions needed), that works in a Euclidean 
domain. 

We then ask some questions related to the results, and mention some 
similar questions that have been answered. In the final section (written 
jointly with P. Moree), some quantitative statements regarding the asymp- 
totic behaviour of various sets of primes satisfying related properties, are 
considered. 
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Part 1 : Divisibility Tests 
1. The two divisibility tests: going forward and backward 

We are all familiar with divisibility tests for a few integers such as 3, 9, and 11. 
To test whether a number is divisible by 3 or 9, we look at the sum of its digits, 
which is equivalent to taking the weighted sum of the digits, the weights all being 
1. For 11, the corresponding test is to examine the alternating sum and difference 
of digits, which means that the weights are 1,-1,1,-1,... . 

The problem of finding a sequence of weights for various divisors has been dealt 
with in [Rh ; in this section, we briefly mention the various tests. First, some 
notation. Any s £ N with digits Sj : < j < m is a polynomial of 10, i.e. 

in 

s = s m s m -i . . . s 1 s = ^2 S J WJ =: *s( k ) 

3=0 

and similarly, one defines s(k), given s,k & N. The coefficients are the digits, i.e. 
< Sj < 9 for all j. To test whether s is divisible by d 6 N, let k = 10 mod d, \k\ 
minimal. Then 

s = s (10) = s (k) mod d 

and thus s is divisible by d if and only if s (k) is, i.e. the sequence of weights is 
simply k , k 1 , k 2 , . . . , where 

k = k d = 10 mod d. (1) 

(In any base B > 1, this test would correspondingly use k d = B mod d.) 

This proves the divisibility tests for 3, 9, and 11 (the last one holds because 
10 = —1 mod 11). It also helps test for powers of 2 and 5 (e.g. a number is di- 
visible by 25 iff its last two digits are 00, 25, 50, or 75). However, if n > 20, then 
k = 10 mod n, and the one with smallest \k\ is 10 itself. So this does not simplify 
the testing of divisibility. 



We now mention the reverse test, also given in |Khj . Suppose d € N is now 
coprime to both 2 and 5. Then 10 has a multiplicative inverse modulo d, and we 
shall take k to be the representative of least absolute value. Thus, 10/c = 1 mod d 
and \k\ is minimal. Then given s = s m s m _i . . . siSq = Y^JLo SjlO 3 , we have 

m m s m \ 

s = j^ sjio j = y, s i k ~ j = wm ( Yl s i km ~ j ) mod d - ( 2 ) 

3=0 3=0 ^j=0 ' 

Since gcd(d, 10) = 1, s is divisible by d if and only if 

in 

7(k) = Y s J km ~ j 

3=0 
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is too. Thus, the sequence of weights (proceeding right to left this time) is simply 
A: , k 1 , k 2 , . . . , where 

k = k d = HT 1 mod d. (3) 

Observe that we can also use equation (J2J) to compute s mod d even when s is not 

divisible by d. We simply reduce 10 m mod d and multiply the result by s(kd) 
mod d. 

For example, A; 29 = 3. Thus, to test 841 for divisibility by 29, we compute 
841(3) = 8 + (4 x 3) + (1 x 3 2 ) = 8 + 12 + 9 = 29, which is divisible by 29. Hence 

29 1 841. Note that this test is better than the standard one, because ^29 = 10 or 
—19. (As a later result on fc-values for arithmetic progressions shows, we might, of 
course, get large values of k for large values of d.) Note that the tests for 3, 9, and 
11 work this way just as well. 



2. The proofs for a general Euclidean domain 
2.1. Preliminaries. 

Definitions. 

(1) Let (E,v) be a (commutative) Euclidean domain with unity. Thus, E is 
a commutative integral domain, with unity and a valuation v : E \ {0} — » 
N :=NU{0}, such that: 

(a) Vs, 7^ d in E, we can carry out the Euclidean algorithm. In other 
words, there exist q, r € E so that s = dq + r, with r = or v(r) < u{d). 

(b) If a, b E E are nonzero, then v(ab) > v(a). 

(By general theory, E is also a PID, and hence a UFD.) 

(2) Define E := E x U {0}, where E x is the group of units in E. (Thus, E = E 
iff i? is a field.) 

(3) Given (s,d) as above, define go = s. Given q n , inductively define g n+ i and 
r n to be any choice of quotient and remainder, when one divides by d. In 
other words, q n = dq n+ \ + r n for all n. 

We say that the Euclidean algorithm terminates for (s,d), if there is a 
sequence of q n 's, that is eventually identically zero (equivalently by Lemma 
n below, some q n = 0). 

(4) We also say that the strong Euclidean algorithm holds if given s £ E and a 
nonzero nonunit d £ E \ Eq, for any choice of quotient q and remainder r 
(in the Euclidean algorithm), we have u(q) < v(s) (resp. v(r) < v{d)), if q 
(resp. r) is nonzero. 
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For example, the strong Euclidean algorithm holds for E = Z and i^f-X"], F a field. 
The latter is easy to see, since the degrees of q, r must be less than that of s, d 
respectively; for the former, suppose s = dq + r. If \dq\ < \s\ then we are done; else 
\dq\ > |s|. But then r / 0, so \q\ = \s - r\/\d\ < (\s\ + \r\)/\d\ < (\s\/\d\) + 1. But 
this is at least (|s|/2) + 1 since d ^ Eq, so \q\ < (|s|/2) + 1. If \s\ > 2 then we are 
done; otherwise d > 2 = s, and we have 2 = d ■ 1 + (2 — d) = d ■ + 2. In both cases, 
u(s) = 2 > u(q). 

We now show two easy lemmas; the first shows what it means to satisfy one of 
the defining axioms of a Euclidean domain, and the second explores properties of 
units, related to the valuation. 

Lemma 1. Given a commutative integral domain E with unit, that satisfies the 
Euclidean algorithm, (as in (l)(a) above) for some v : E — > No, the following are 
equivalent: 

(1) The Euclidean algorithm gives unique quotient and remainder, when applied 
to (de, d) for any e G E and any d ^ 0. 

(2) The Euclidean algorithm gives unique quotient and remainder, when applied 
to (0, d) for any d ^ 0. 

(3) v(ab) > u(a)\fa,b^0. 

Proof. Clearly, (1) implies (2). Given (2), suppose v(ab) < u(a) for some a, b. Then 
if we divide by a, we have two solutions to the Euclidean algorithm, namely 

= • a + = (-6) • a + (ab) 

which is a contradiction. Finally, given (3), suppose we have a solution to de = 
d-q + r. Thus r = d(e — q), so that if e ^ q, then i/(r) > u(d) > v{r), a contradiction. 
Hence e = q and r = 0. □ 



Lemma 2. The following are equivalent for 7^ B € E: 

(1) B G E x (i.e. B is a unit). 

(2) v(Be) = v{e) for all nonzero e G E. 

(3) v{B) = v(V). 

(4) v{B) < v(l). 

(5) The set {0 / o G B : u{a) < v(B)} is empty. 

Proof. Firstly, if B is a unit, then u(e) < v(Be) < v{BB~ 1 e) = v{e), so they are all 
equal and (1) => (2). Setting e = 1, (2) =J> (3) => (4). Next, (4) implies (5) because 
v{a) = v(a ■ 1) > ^(1) > v(B) for all nonzero a, so that the desired set of elements is 
empty. Finally, assume (5). Then 1 = qB + r by the Euclidean algorithm, and r = 
or u(r) < v(B). Given our assumption (5), we conclude that r = and 1 = qB, so 
that B is a unit. □ 

We now adapt the results of the previous section to (E, v). Henceforth we assume 
that E is not a field. Fix a base B ^ Eq. From Lemma |U above, this means that 
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we can talk of "digits" for numbers. Fix also a divisor d G E coprime to B (i.e. 
g.c.d.(d,5) = l). 

Lemma 3. 

(1) If the strong Euclidean algorithm holds, then it terminates for all (s,B). 
In other words, every s G E can be written as s = s m s m -i . . . s\So = 
Y^jLo s jB 3 , where Sj = or v{sj) < v{B) for all i. 

(2) If(B,d) = I, then 3\ B" 1 mod d. 

Proof. The first part follows from carrying out the Euclidean algorithm repeatedly, 
i.e. repeatedly dividing by B. Note that at each stage, the quotient q n +i has a 
strictly smaller v- value than the previous quotient q n , whence it must eventually be 
less than v{B). We then define q n +2 = and r n+ \ = q n +i- 

The second part holds in any PID: since B and d are coprime we can find oq, /3q G 
E such that a$B + fad = 1. Thus ct$B = 1 mod d, and olq is an inverse modulo d. 

If there are two such, namely -B-f and -B^~ , then (£>j~ — B2 )B = 1 — 1 = 
mod d, so it equals Id for some I G E. Multiplying by oq, and noting that a^B = 
1 — Pod, we get 

(B^ 1 - B^ 1 ) = [(B^ 1 - B^)/3 + la ]d 
which means that d\LHS, and B^ 1 is unique modulo d, as claimed. □ 



2.2. The proof of the reverse test. 

We again test whether or not s is divisible by d in E. Clearly, we have two dif- 
ferent methods, as above. The first (or standard) method is as above, and we focus 
on the second one now. For notation we fix B > 1, our base, and d, our divisor, 
that is coprime to B. Given a dividend s G E, let sq be the unit's digit and (3 be 
the "rest", i.e. s = f3B + so, with sq = or u(sq) < v{B). 

Theorem 1. There is a unique k G E/dE, that satisfies the following condition for 
all s G E : 

s = j3sq = mod d iff (3 + kso = mod d 
Moreover, 

k = k d := k d = B^ 1 mod d (4) 

Proof. Existence is easy: define k d = B^ 1 mod d. Now we see that f3 + k^SQ = 
kd{Bf3 + so) — kds mod d. Since B and d were coprime, hence so are kd and d, so 
s = Bj3 + so = mod d iff f3 + k^so = mod d, as desired. 

For uniqueness, suppose any such k exists. We then keep kd as above. Now take 
any multiple s of d, so that so = 1, i.e. s = 1 mod B. Since d and B are coprime, 
such a multiple clearly exists (e.g. from the proof of Lemma El above, Pod = 1 
mod B). Then by our condition, /? + kso = P + k = mod d. But since s = 
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mod d, hence k^s = mod d, whence we get k^s = kd(B/3 + 1) = f3 + kj = 
mod d. 

Thus (3 + kd = f3 + k = Q mod d, whence k = kd mod d, as required. □ 

Notation : Henceforth, we only consider dividends of the form s = s m . . . so, where 
the Sj's are the digits of s, when written in base B. Let < n\ < ■ ■ ■ < n\ < m for 
some I. Then define Oi nit ,., i n ha {k ( i) to be the sum 

c a Xt" 1 ""' 1 - 1 , , o A_k ni ~ n '- 2 ^ 1 S -I \-k n '? 1 In 

°m ■ ■ ■ °ni i ">(} a nj — 1 ■ • • °ni — i ' "'d a ni—i— 1 ■ ■ ■ °ni_2 > < "'d a "-i _ 1 • • • °0 



where s a . . . st, is the number with digits s a , . . . , Sb respectively (when written from 
left to right, in base B). 

Corollary 1. cx ni) ... >ni ,s(kd) = B~ n 's mod d. 
Proof. We compute, for general i, 

rti+l — l Mj+i— 1 Tlj+1 — 1 

K d Sn l+1 -1 ■ ■ -Sm — K d / 4 SjO = ^ S 3 K d — /_> 3d 

j=m j=n t j=m 

Therefore we have 

m m m 

a ni ,..,nM - E s i k T j = k 7 E *j V - B~ ni E S > BJ = B ~ nis 

3=0 j=0 j=0 

(where everything is modulo d). Hence we are done. D 

Remarks: Thus, the number kd (as above) also satisfies (for any s) : 

(1) For all such tuples (ni, . . . ,rii), s = mod d iff u ni ,...,rn,s(kd) = mod d. 

(2) Hence for all i, s = mod d iff ai >s (kd) = s m . . . Sj + k d Si-\ . . . Sq = 
mod d. 

(3) Setting m = i Vi = 1, . . . , m — 1 in (1), we get that s = mod d 

iff ai t 2,...,m-i,s{kd) = a s (k d ) = mod d. 

In particular, the "reverse" divisibility tests (in Theorem^and Section^) mentioned 
above, do hold. 

Example: There is a well-known divisibility test for 7 that involves splitting up 
numbers into groups of three digits, and then taking the alternate sum and difference 
of the numbers. For instance, to check 142857142 for divisibility by 7, we compute: 
142 - 857 + 142 = 284 - 857 = -573 = -13 = 1 mod 7, whence 7 \ 142857142. The 
reason this works is the corollary above (and the remarks following it), if we note 
that 1000 = -1 mod 7, so that 10 6 = 1 mod 7. 



DIVISIBILITY TESTS AND RECURRING DECIMALS IN EUCLIDEAN DOMAINS 7 

We now show that the /c-values of terms in an arithmetic progression (more pre- 
cisely, having the same "unit's digit"), themselves form an arithmetic progression. 
This is useful to calculate /c-values for general d, if we know the /c- value of only the 
unit's digit of d. (For example, knowing kg, we can easily calculate k\g, /C29, etc.). 

Theorem 2. Take a E E (a,B coprime), and m E E. Let d = Bm + a. Then 

kd = ka+mB = k a + ml mod d, where I = —(Bk a — 1) (5) 

a 

Proof. We compute: B(k a + ml) = Bk a + Bml = {l + al) + Bml = l + l(a + mB) = 1 
mod (a + mB) = 1 mod d. Hence by definition of kd, we are done. □ 



Lemma 4. Let d\(B f — a) where a E E x and t E No- Then kd = B t (aB) x mod d. 

Proof. B{B t {aB)~ 1 ) = eT 1 ^*) = a~ 1 (B t - a) + 1 = 1 mod d = Bk d mod d. 
Hence we cancel B (modulo d) to get the result. □ 

This (last) result tells us, in particular, the divisibility tests for 3, 9, and 11. We 
now give two examples in the Euclidean domain E = F[X] (where F is a field.) 

Example 1: Let us work with B = X, d = X 2 + 1. We want to check whether or 
not s = X 4 + X 2 is divisible by d. Clearly, s = 10100 in base X, and to apply our 
test, we note that (— X)d = —X 2 = 1 mod d. Hence k = —X = —10. So we get 
a s (k) = 1 + (— 10) 2 = 101 and this is indeed divisible by 101 = d. Therefore, so is 
the original s. 

Example 2: We now show how the Factor Theorem for polynomials is an example 
of the reverse test. We still work with B = X, and d = X — c for some c E F x . We 
want to check if a general element 

m 

s = s m s m -i ...siSq = ^2 s j XJ 

3=0 

is divisible by d or not. From Lemma 0] above, kd = c~ x mod d. Thus to check 
whether or not d\s is the same as checking whether or not (X — c) divides a s (kd) = 
YlT=o s j( c ~ 1 ) m ~'' ^ F- Since this is a scalar, it should be zero, and multiplying by 
c m should still give zero. Thus, 

m 

(X - c)\s(X) iff J2 S J^ = ° iff s ( c ) = ° 

3=0 

and this is just the Factor Theorem for polynomials. Moreover, Theorem [21 above 
says that the "fc-value" function (taking d m = mX — c to k(d m )), is a linear function 
in m E F. 
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2.3. Testing for divisibility by general n E E. 

We now check whether or not a given number s is divisible by n, where n E E 
need not be coprime to B, the base in which both s and n are written. To do 
this, we take n = dp^p^ 2 ■ ■ -Vy^ (by unique factorization in the UFD E), where 
d is coprime to B, and pi,P2, ■ ■ ■ ,Ph are some or all of the primes that divide 
B (on EN for i = 1,2, ...,h). 

Divisibility by d is tested by the method given above, and for each pi, we perform 
the analogue (in base B) of the "n-right-digits test" : 

If Pi is the highest power of pi dividing B, then to test for p?*, we check whether 
or not the number formed by the |~aij//3j] rightmost digits of s (in that order and in 
base B) is divisible by p^\ (Here, \x] is the least integer greater than or equal to 
x, e.g. if /3i\cti, then \ai/Pi\ = ati/fti.) This is because the place value of each of 
the other digits is > B' 01 ''^'', and this is already divisible by pr 1 ' 011 ' * ' ; which is a 
multiple of pf* . 

For example, take E = Z. In base 10, to check for divisibility by 125, we only 
need to check the three rightmost digits - they have to be 000, 125, 250, ..., 750, 
or 875. As another example, consider B = 12. Since 12 = 3 x 2 2 , hence to check 
divisibility by 8 = 2 3 , the [3/2] = [1.5] = 2 rightmost digits have to be checked. 



Part 2 : Recurring Decimal Representations 

3. Decimal representation of fractions - a few observations 

We know that every rational number can be represented as a number having 
a terminating or a recurring set of digits in its decimal representation. We ob- 
serve certain interesting properties of the digits in the recurring portion, e.g. j = 
0.142857 142857 . . . , and upon multiplying 142857 by any number from 1 to 6, we 
get the same six digits in the same cyclic order (142857, 285714, 428571, 571428, 
714285, 857142). On multiplying it by 7, we get 999999, so that in the dec- 
imal representation, it becomes 0.9999- •• = 1. Another similar case is ™ = 
0.05882352 94117647 0588..., and if 0588 ... 7647 (including the zero) is multi- 
plied by any number from 1 to 16, then the same sixteen digits are obtained in the 
same cyclic order; if it is multiplied by 17 then the result, of course, is 9999 . . . 9999 
(sixteen 9's). This is further addressed in the comments following Lemma ITUl b elow . 



We also observe the following phenomenon, which we quote from |Kvj : We know 
that 1/7 = 0.142857. If we square 142857, we get 20408122449. Take the first six 
digits (from the right) and add them to the number formed by the rest (we take six 
digits because 142857 has six digits). We then get 20408+122449 = 142857. 

Similarly, if we take 1/27 = 0.037 and square 37 then we get 1369, and adding the 
rest to the three rightmost digits yields 1+369 = 370, again a multiple of the original 
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number. We shall also show this phenomenon in general, in Lemma ITUl below. 

We observe another interesting property in both cases - and more generally, in the 
case of every prime whose reciprocal has an even number of digits in its recurring 
decimal portion. Namely, if the second half of numbers is kept below the first half, 
then the sum of every pair of corresponding digits is 9 (e.g. in 142 857, 1+8 = 4+5 
= 2 + 7 = 9). We give further examples below, and explain this in Lemma El and 
Theorem |3] below. 

First of all, how do we get recurring decimals ? For an odd prime p / 5, we 
look at the order b of 10 in (Z/pZ) x . We have W b = 1 mod p. Thus p divides 
10 6 - 1 = 999 ... 9 (b nines). Let 10 6 - 1 be denoted by %, and let 2fc = Tp , Thus, 

we have jjj = ^ (0.999 . . . ) = |(0.9fe9fe . . . ) = 0.r p r p Then the number of digits 

in the recurring part of the reciprocal of p gives the number of terms (i.e. b) in the 
chain. Related facts are shown in Lemma [7| below. 

Given p, we have various chains of recurring decimals - for instance, we have the 
(unique) chain 142857 = 571428 etc. for 1/7. However, for some numbers (e.g. 13), 
we might have more than one chain. Let c p be the number of distinct chains. Then 
b p Cp =p-l. 

Now consider the situation in base B (for any B > 1). In what follows, the 
phrase "decimal expansion" also refers to the case of a general base B. In general, 
for composite d, the recurring parts of t/d might have recurring chains of different 
lengths as well (e.g. for d = 21, in base 10, the chains for 7/21 and 14/21 have length 
1 each, and are different from the 6-digit chains for the other residues). However, if 
we denote the different lengths of all possible chains by bi , and the number of chains 
of a given length (= bi) by Cj (for i = 1, 2, . . . ), then we have 

For any d (d € N, d coprime to B), N^ b{Ci = d (6) 

i 

We also remark that the number of chains c p is also known as the residual index, 
and it clearly equals [(Z/pZ) x : (J3)]. 

Note that in base 10, we take d/d = 1 = 0.999 • • • = 0.9, and this makes the above 
sum one more than the d — 1, that we obtained above, in the case when d was a 
prime. Similarly for base B. 

Of course, we can find many composite numbers d which have only one value of 
b (= bd, say), e.g. 691 = 6. This happens, for instance, when d is the product of two 
or more primes, all distinct, and all having the same 6-value. This is addressed in 
Proposition ^ below. 

Proof of Equation (jfJJ). Look at {t/d : 1 < t < d}, and look at their "dec"imal 
representations. Clearly, there are no repetitions among the various decimals, for the 
numbers in the above set are all distinct (modulo d). Further, the numbers in the set 
whose recurring chains coincide with the chain t/d = 0.a\ . . . a^ are precisely t/d, (Bt 
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mod d)/d, (B 2 t mod d)/d, . . . , (B b ~ 1 t mod d)/d. Hence, given any chain, every 
cyclic permutation of it corresponds to t/d for some d. 

This gives a bijection between all cyclic permutations of a fixed chain (6j of them), 
and bi of the residue classes modulo d. Hence we get equation © above, because 
there are d residue classes modulo d (where we also include 1 = d/d = 0.9 etc.). □ 



4. Properties of recurring decimals - proofs for a Euclidean domain 

We now work on fractions and recurring "dec" finals in base 5, in a Euclidean 
domain (E,v). Though we cannot talk about convergence here, we can still talk 
about repeating decimals representing a fraction, and such things. 

Observe that if ^ d is not a unit, then z/(l) < v{d), so 1/d should have an 
expansion of the form 0.a±a2 • • • = ^ a>iB~ l € k(E), the quotient field of E. (Actu- 
ally, we work in i?[[i3 _1 ]] = E[[X]]/(BX — 1).) We now find out what the Oj's are. 

Remarks : 

(1) We do not worry about convergence issues here, and only associate to a given 
fraction, a sequence a-i of numbers, without worrying whether these numbers 
are actually digits or not (i.e. whether u{ai) > or < v{B)). Note that if 
E = Z, then all sequences are actually decimal expansions that converge etc. 

(2) If d ^ and d \ a in E, then a/d = Q.a\a2 ... iff v{a) < v{d). In other words, 
there exists a representation of a/d as 0.aia2 ... iff the quotient (when we 
divide a by d in E) is zero, iffa = 0-d + r = r satisfies the Euclidean 
algorithm, namely, that v(a) = u(r) < v{d). 

Lemma 5. Suppose ^ d is not a unit, and we define h = 0.a±a2 ■ ■ ■ ■ Then 

a l =[B i /d\-B[B i - 1 /d\. 

Here, by [a/b\ , we mean the quotient when we divide a by b. 

Warning : We need to choose and fix the quotients [B l /d\, since these are not, in 
general, unique ! Moreover, convergence cannot be discussed, and the cjj's need not 
be digits (i.e. we do not know if we have v{ai) < v{B) or not). 

Proof. Suppose we write 1/d = 0.a\a2 .... Clearly, the quotient (i.e. the number to 
the left of the decimal point) is zero because of the Euclidean algorithm (and since 
i/(l) < v{d)). Next, a\ is clearly of the desired form, because of the same reason. 

We now show that the a,'s are as claimed, by induction on i. The base case was 
done above. To show the claim for «j, multiply the equation above by B l . Thus, 
B l /d = a\a2 ■ ■ ■ ai.ai + \ai + 2 . . . , so taking the quotient of both sides yields 



YB l /d\ = a l a2...a l = YB i 
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This is because B l = d{a\ . . . Oj) + d(0.a i+ i . . . a n a\ . . . a^) = dq + r, say (where 
r = d(0.aj + i . . . a n a\ ... a,) 6 -E). By the above assumption (since (B,d) = 1), we 
get that v(r) < v(d), so q and r are clearly the quotient and remainder, as desired. 

Assume by induction that we know the results for ax,... , a%—i. Then 

i-l 
a { = [B l /d] -J2 Bi ~ J (l BJ /d\ ~ B[B^ x /d\) 

i=i 
and the latter is a telescoping sum, so we get 

at = [B l /d\ -B[B l - l /d\ +B\B°/d\ 

But the last term is zero, because u(B°) < v(d). Hence we are done. □ 



Lemma 6. The recurring decimal O.ai . . . a n denotes the fraction ax ■ ■ ■ a n /(B n — 1) . 
If a\ . . .a n = B n — 1, then the decimal equals 1. 

Proof. Suppose e = O.ai . . . a n . Then B n e = ax ■ ■ ■ a n .ax ■ ■ ■ a n . Subtracting from 
this the definition of e, we are done. If, moreover, ax ■ ■ -a n = B n — 1, then (B n — 
l)(e — 1) = 0. Since B is not a unit, B n ^ 1; since E is an integral domain, we 
conclude that e = 1. □ 

Next, to talk about periodicity of any and every a/d (for fixed d coprime to B) 
in such a decimal expansion, we need the following 
Standing Assumption : B has finite order in {E/dE) x . (In particular, d ^ E x .) 

Denote this order by e. Thus e = Od(B) = Od{B~ l mod d) = o^kd). We now get 

1 B e - 1 . 1 1 N , B e - 1 

— = ( 1 k — h • • • ) = O.rjr^ . . . , where rj = 

d d y B e B 2e ' ' d 

and this is how we get recurring decimal expansions. Moreover, from the proof of 
Lemma above, r^ = (B e — V)/d = [B e /d\ (because u(l) < v{d)) = ax ■ ■ ■ a e , and 
hence r^ has e digits. Thus, the correct way to look at this, in a Euclidean domain, 
is to look at the order of B in {E/dE) x , instead of the recurring decimal expansion 
or the number of digits therein. 

Lemma 7. If t is coprime to d, and B has finite order e in {E/dE) x , then the 
sequence of a, 's associated to t/d is recurring with period e. 

Proof. The length of the sequence is the smallest e' S N such that tB e = t mod d. 
Since t is coprime to d, we get e' = e by definition. □ 

We now turn to repeating decimals with even period. We have 

Lemma 8. B has (finite, and) even period, modulo a prime p \ IB, iff ' p\(B l + 1) 
for some I € N. 



12 APOORVA KHARE 

Proof. Suppose p\(B l + 1) for some positive I; we assume, moreover, that I is the 
least such. Thus p\(B 21 — 1), whence the order o p (B) is finite and divides 21. Since 
p \ 2, hence p \ B — 1, whence o p (B) \ I. Thus o p (B) = 2m for some m|l. 

Conversely, if o p (B) = 2m, then p\(B m + l)(i3 rri — 1). By definition of order, 
p \ [B m — 1), whence p\ (B m + 1), and we are done. (Note, moreover, that continuing 
the above proof of the first part, we find that I = mby choice of I.) □ 



Lemma 9. Suppose d\{B l + 1) for some I € N. Then for all a coprime to d, the 
chains starting from the ith and the (I + i)th "digits" in the recurring part of a/d, 
add up to a recurring chain that represents 1 or 0. 

Proof. Since B l = — 1 mod d, hence the chain corresponding to adding up the frac- 
tions (aB l mod d)/d and (aB l+% mod d)/d, yields the fraction {(aB % + a(—l)B l ) 
mod d)/d = d/d = 1, as claimed. □ 

Examples: 



(1) Suppose E = Z and B = 10, so that v{B) > u(B-l). Then 1/7 = 0.142 857, 
which satisfies the above assumptions, since the chains add up to 0.999999 = 
0.9 = 1. __ 

Similarly, 1/11 = 0.09. Here's another way to write this. Observe that 
9 = 10-0 + 9 = 10-1 + (-1), with i/(-l) = | - 1| < i/(10). Hence we can 
also write 1/11 = 0.1(— 1), and then 1 + (—1) = 0. This corresponds to the 
recurring chain 1.0 for 1. 

(2) As another example, consider B = F[X], with char F/2. Suppose B = X 
and d = X 2 + 1 (so that v(B) = v(B — 1)). Clearly, B has order 4 modulo 
d, since (i|(X 4 — 1). So 1/d should have four digits. But we also know what 
these are : r d = (B e - l)/ d = (X 4 - 1)/{X 2 + 1) = X 2 - 1 = 010(-1) in 
base B = X. Thus 1/d = 0.010(— 1) and once again, every pair of alternate 
digits adds up to zero. 

We thus see that the case E = Z and 1 = 0.9 is somewhat of an "exception" ! This 
is made more precise now. 



Theorem 3. Given a Euclidean domain (E, v) that is not a field, the following are 
equivalent: 

(1) u{B) = u{B-l) VB (£E . 

(2) v(B) = v(B - c) | VS i Eq, c e E . 

(3) The Euclidean algorithm yields a unique quotient and remainder when we 
divide any B n — c by B, for any n £ N, B ^ Eq, c £ Eq. 

(4) There is no B ^ Eq, so that 1 can be expressed (in base B) as a recurring 
decimal O.ai . . . a n for some nGN, with v{aj) < v{B) or a« = 0. 

Proof. We prove a series of cyclic implications. 

(1) => (2): 



DIVISIBILITY TESTS AND RECURRING DECIMALS IN EUCLIDEAN DOMAINS 13 

This is clear for c = 0, and for c £ E x , we have (assuming (1) and using part (2) of 
Lemma |2j) 

V (B - c) = v{c~ x {B - c)) = v{c~ x B - 1) = v(c~ l B) = v(B) 

(2) => (3): 

Suppose -B n — c = (?£> + (i. Then B\(c + d), whence (c + d) is not a unit. If d ^ — c 
then c + <i is a nonzero nonunit (whence <i 7^ 0). Hence z^(c + d) = v{d) by (2). 
Moreover, B\(c + d), so f(-B) < u(c + d) = v(d) < v(B), the last inequality a conse- 
quence of the Euclidean algorithm. This gives a contradiction, whence d = — c and 
hence q = B n , as desired. 

(3)=>(4): 

Suppose we can write 1 as a recurring "decimal" . Then by Lemma B n — 1 = 
a\ . . . a n = a n +B-a\ . . . a n -i- Assuming (3), one notes that a n = —1, since —1 & E . 
Thus a\ . . .On-i = B n . Once again, by (3) we get that a n _i = 0. Proceeding 
inductively, we get that a, = VI < i < n, and finally, we are left with a\ = B, a 
contradiction to v{a\) < v(B). Hence 1 cannot be written in recurring form. 

(4)=>(1): 

We prove the contrapositive. Suppose we have a B £ Eq so that v(B) 7^ v{B — 1). 
There are two cases. If v{B) > v{B — 1), then (using Lemma EJ) we have the 
expansion 1 = 0.(B — 1) in base B. On the other hand, if v{B) < v(B — 1), then 
the expansion e = O.B(-B) in base {B — 1), again equals 1 (by LemmaEJ). □ 



Lemma 10. Suppose B e — 1 = d ■ r^, where r^ is the recurring part. Fix k, I € N 
and consider krd- Then the number formed by the first le "digits" (or ai's), added 
to the number formed by the rest, yields a number that is divisible by r^. 

Proof. The number in question is obtained as follows: suppose the "rest of the digits" 
form a number s'. (For instance, if E = Z, then we could write s' = [krd/B le \ , 
where this denotes the greatest integer in the quotient.) Then the total number 
obtained is 

s' + {kr d - B le s') = kr d - (B le - l)s' 
and this is divisible by r d , the quotient being k — ds'(l + B e + • • • + B le ~ e ). □ 

We now show some of the facts from the previous sections. First of all, observe 
that if t/d = O.ai . . . at . . . , then 0.aj + iaj + 2 • • • comes from the expansion of B l t/d, 
or from {B l t mod d)/d = a/d, say. Thus if r d = (B e — l)/d, then the recurring part 
of a/d would correspond to ar d . This explains why recurring decimal expansions 
keep getting permuted cyclically, upon being multiplied by various numbers. 

Next, we define D d , for any d satisfying the above assumption, to be the order e 
of B in {E/dE) x , or equivalently, the "number of terms" of the sequence (at), in r d 
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above. 

Suppose d = \Yl = iP^ i as above, where pi are primes in E (which is a unique 
factorization domain). If d'\d, then 

-> dE -> d'£ (--» £7) 

and hence 

£7/dE -» E/d'E -» 
(The kernel is, of course, d'E/dE, so we get the short exact sequence 

-» d'^/dS -> £/d£ -» E/d'E -» 

for every choice of d' .) In particular, B has finite order in every (E/d'E) x (because 
if B e = 1 mod d, then B e = 1 mod d'); in particular, in {E/piE) x as well. Thus 
we can talk of D Pi as well. However, this does not guarantee that the standing 
assumption holds for p" for all n. So we need to make a stronger assumption. 

However, let us remark that once B has finite order in (E/dE) x , every t/d has 
a recurring decimal expansion. For, we write t/d = if jd! in lowest terms, and then 
from the above remarks, B has finite order in (E/d'E) x , and we are done. 

Proposition 1. Suppose d = p\ . . .pt for distinct primes pi, all coprime to B. Let 
bi (£ N) denote the order of B in (E/piE) x . 

(1) Then B has finite order in {E/dE) x , and this order, denoted by bd, say, is 
the 1. cm. of the bi 's. 

(2) If bi = bd Vi, then B has order bd in {E/d'E) x , for every nonunit d'\d. 

Proof. 

(1) Clearly, B lcm — 1 = mod pi Vi, where 1cm denotes the l.c.m. of the ftj's. 
Since the p^s are mutually coprime, and E is a UFD, hence YliPi also divides 
B — 1. Thus B has finite order in (E/dE) x , and moreover, this order has 
to divide the l.c.m. But the l.c.m. has to divide the order as well, so we are 
done. 

(2) Suppose bi = bd Vi. Now, if d \ a, then a/d is of the form r/s, where r and 
s are coprime, and s is a product of some of the p^s. Clearly, r/s and 1/s 
have the same number of recurring digits, because (r, s) = 1. But s ^ 1 is a 
product of some p^s, and each 1/pi has bd recurring digits. Hence so does 
1/s, and so does r/s = a/d. 

□ 

5. The period of recurrence of a fraction in a Euclidean domain 

We have the following 

Lemma 11. If B has finite order e modulo p, then the following are equivalent: 
(1) B has finite order modulo p n for any n. 
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(2) There exists a unique prime integer p' G Z (which is prime in 7L, but may 
not be prime in E), so that p\p' ' . 

(3) There exists n G Z so that p\n. 

Remark: By Z we mean the image of Z in E, for we do not know if E has charac- 
teristic zero or not. 



Proof. Suppose B has order e modulo p. Thus B e = 1+ap for some a G E. Suppose 
a = p*(3, where p \ /?. Now, B en = 1 + nap mod pf +2 (ignoring higher order terms). 
Thus, B has finite order modulo p* +2 iff there is some n such that pf +2 \nap. But 
from our assumptions on /, this means p\n. Thus (1) => (3). 

Conversely, assume that p\n. Now observe that if s = 1 mod p m for any s G 
E, m G N, then s = 1 + /3p m for some (3 G -E 1 , whence 

s ra = (1 + /?p m ) n = 1 + [3np m mod p 2m = 1 + [3np m mod p m+1 = 1 mod p m+1 

since p\n and m G N. In particular, since B e = 1 mod p* +1 = 1 mod p, hence 
B enm ^ = 1 mod p m for all m G N. Hence (3) => (1). 

That (2) and (3) are equivalent is easy : in one direction, choose n = p'. In 
the other direction, decompose n into prime factors (as integers), and note that 
in E, p must divide some integer prime factor of n. If p divided two such, then 
it would divide any integer linear combination of these prime numbers, whence p 
would divide 1, and hence be a unit. This is not possible. □ 

Thus, we now make the following 

Standing Assumptions : 

(1) (E, v) is a Euclidean domain, not a field. Fix B $_ Eq. 

(2) For each prime p G E, p\ B, (a) B has finite order modulo p, and (b) there 
exists a unique prime integer p' G Z so that p\p' (and p' ^ 0) in E. 

As an easy consequence, we have 

Lemma 12. E (or its quotient field) has characteristic 0. 

Proof. E is not a field, so there exists a nonzero nonunit in E. By unique factoriza- 
tion, there exists a nonzero prime p G E, so by assumption there exists a nonzero 
prime integer p' that is a multiple of p. Now, p' ^ E$ since p £ Eq. This would give 
a contradiction if E had positive characteristic, so we are done. □ 



Lemma 13. If (B,d) = 1, then B has finite order modulo d. Further, there is a 
unique n G N such that for each meZ, d\m in E iff n\m in Z. 

Proof. We argue by induction on the number of distinct prime factors of d (ignoring 
multiplicities). If this number is 1, then d is a prime power, and we are done by the 
previous lemma. 

Suppose we now have a general d = p r d' for some p prime in E, coprime to d' , 
and suppose we know the result for d! '. Then B* = 1 + ad' for some a, /. If B e = 1 
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mod p r ', then B e $ = 1 modulo both d! and p r , which are coprime. Hence by unique 
factorization, B e -> = 1 mod d, and we are done. 

Finally, if d = UiP?, and ft|pj e Z, then dWl^p'^ € Z. Thus the set S of 
integers divisible by d, properly contains the element (since E has characteristic 
zero, from above). Moreover, S is clearly an ideal in Z, hence is principal and 
generated by some n£ff. This is the required n (and d\n), and we are done. □ 



We now calculate D^ = the number of digits in the recurring decimal part of §, 
where d is coprime to both x and to B, and d = p^p^ 2 ■ ■ -pV 1 (the p^s here are 
distinct primes coprime to B, and a, £ N). 

For the rest of this section, fix d = Yl^pf 4 coprime to B, and fix p, a prime factor 
of d. Define qi{p) = q% = D p i = order of B in (E/p l E) x , and define / to be the 
highest power of p that divides the prime integer p' (as in the Standing Assumptions 
above). Also let v p {m) be the largest power of p that divides m in the UFD E. 

We first prove the following generalization of Ri, P.1.2 (iv), pp.11] (special cases 
of the result in [EQj were known to Euler) : 

Proposition 2. Given x ^ y in E, define a(m) = - ^ for m > 0. Now assume 
thatp\(x—y) andp\y. Thenp\a(m) •£> p'\m. If, moreover, v p (x—y) > [f /(p'—l)\ € 
No, then for all m6N, we have 

u p (a(m)) = Vpl = v p {m) 

\ x-y J 

In particular, (x p — y p )/{x — y) is divisible by p* but not by p* +l . 

Proof. Firstly, it is easy to check that (x—y)\(a(m)—my m ~ 1 ) for all m. Now suppose 
p\(x — y). Then a(m) = my m ~ l mod (x — y), whence the same holds modulo p. 
Since p, y are coprime, we conclude that p\a(m) 44> p\m <=? p'\m, by Lemma 1131 
above. 

Now suppose that u p (x — y) > f/(p' — 1). Setting x — y = t, we get that 

a (p') = {(y + t) p ~ y p )/t = Yli=i Vi)y p ~*t ■ Over here, every term except 
the first and last one, is divisible by p' ■ p, since there is a binomial coefficient and 
a power of t in each term. Thus, modulo p' +l , we have a(p') = p'y p _1 + t p . 
The last term is divisible by p n ^ p _1 ', where n = v p {x — y) > f/(p f — 1); hence it 
is divisible by p-* +1 . The first term is only divisible by p*, since p,y are coprime. 
Thus the last line of the result is proved. 

To prove the rest of the result, suppose m = up' s , where p' \ u. Now define 
Xt = xW,Yi = y ip,) \ and note that 



x m -y m 
x-y 


_ X s u - Y s u 
X s — Y s 


s-l „ 
1=0 


-Y i+1 


x s u - 
x s - 


V u 

1 s 

-Y a 


s-l Y p' 
1=0 


-Yi p ' 
-Y t 
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where the product on the right telescopes. Since (x— y)|(JQ— 1^), hence u p (Xi— Yj) > 
f/(p' — 1) for all i. Hence by the previous paragraph, each term in the product is 
exactly divisible by p* ; moreover, the remaining term is not divisible by p by the 
first part of this result, since p\u. Hence v p { x x ~jL ) = fa = v p (m). □ 



Theorem 4. For all n £ N, q n +i = Qn or p'q n - Now let g be the least natural 
number so that (a) g > f /{p 1 — 1), and (b) q g+ \ = p'q g - Then for all n > 0, we have 

q g+n = (p') {n/f] q g 

Proof. First, the general case. We assume q n +i ^ qn- Now, B Qn+1 = 1 mod p n+1 = 
1 rnodp™, whence q n \q n +i- So let q n +i = mq n . We now use the first part of 
Proposition |2 setting x = B qn and y = 1. Thus, the order of B qn in (E/p n+1 E) x 
is larger than 1 and divides the prime integer p', whence it is p'. 

Next, let g be as above, and let n G N. By the previous paragraph, q g + n = (p') s q g 
for some s. We now apply the second part of Proposition [21 this time setting 
x = B q v,y = 1. Thus 



X (PT _ yip'Y 



M(p') s ) = f* 



x-y 
and the result follows. □ 

Remarks: 

(1) To conclude this discussion, if d = p^p^ 2 ■ ■ -P^ is coprime to B, then to 
find Dd, we find Qi = q ai (Pi) (i = 1,2,... ,h). Since pi,P2, ■ ■ ■ ,Ph all divide 
B Dd — 1, hence Q%, Qi, ■ ■ ■ , Qh all divide D^. Moreover, D^ is the least such 
positive number. Hence 

D d = lcm(Q 1 ,Q 2 ,...,Qh) (7) 

(2) For E = Z, note that j/ = p is itself prime, so that / = 1, and all our 
other standing assumptions are satisfied, for any B > 1. Moreover, Od(B) 
now does denote the number of digits in the recurring part of the decimal ex- 
pansion (in base B) of 1/d, so we can compute this using the above theorem. 

(3) In the decimal system (B = 10), note that for the first few primes p > 5, we 
have g = 1, though 53 = 2 for p = 3. Thus, D3 = Dg = 1, D27 = 3, D^i = 9, 
and so on, while Dj = 6, -D49 = 6 x 7 = 42 etc. 

One can ask if this phenomenon actually holds for all primes p > 5. In 
other words, are the orders of 10 in (Z/pZ) x and (Z/p 2 Z) x unequal for all 
primes p > 5 ? 

The answer is no. The prime p = 487 satisfies q\ = c/2 = P — 1 = 486 (and 

93 =P{P~ !))• 

However, is there a base B for which this property does hold (for E = Z) ? 
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6. Concluding remarks and questions 

We conclude this part with a few questions: all the properties asked below are 
known to hold or not to hold for E = Z or i^f-X"], and we want to ask whether they 
hold in general, or if we can characterise all Euclidean domains with that property. 

(1) Is there always a "good" quotient and remainder ? In other words, can we 
always find a valuation u, so that the following holds ? 

Given s,0 / d £ fi, we can find q, r by the Euclidean algorithm, such 
that s = dq + r, r = or v(r) < v(d), and with the additional property that 
v{dq) < v{s). 

For example, in E = Z, suppose s = —11, d = 3. Then s = (— 4)d + 1 = 
(— 3)d+ (—2), and we take the good (q, r) to be (—3,-2). 

(2) Does there exist a (sub) multiplicative valuation ? Namely, a valuation v 
satisfying : if u(a) < v(b) and v(c) < v{d), then u(ac) < v(bd). 

(3) Does there always exist a valuation satisfying the triangle inequality ? 

(4) Is there a characterization of all Euclidean domains in which every prime 
divides a prime integer (i.e. p\p' in E, where p' is prime in Z, and nonzero 
in E) ? 

For example, it is true in E = Z and in Z[z] (because a + bi divides 
a 2 + b 2 G Z, so every prime divides an integer, and by Lemma fTTI above, we 
are done). But this property does not hold in E = F[X] (F a field), since the 
polynomial X — 1 is prime, but does not divide any integer (or field element). 

We have a partial answer for this last question (note that we already saw above 
that E has to have characteristic zero here). 

Proposition 3. Suppose (E, v) is a Euclidean domain but not a field. If the strong 
Euclidean algorithm holds, and every prime p £ E divides a prime integer p'GZc 
E, then E is finitely generated as an i?-module over the subring R = Z[i? x ]. 

Proof. We produce a surjection : R[X] — > E, with a monic polynomial in the kernel. 
Note that Eq = E x U{0} C R. We first claim that there is a (set-theoretic) surjection 
: E [X] -» E. 

Since E is not a field, there exist nonzero nonunits. Pick one with the least 
valuation, say B E E, such that f{u) < v{B) < v{y) for all nonzero nonunits y and 
all units u in E. By Lemma |3 the evaluation map : X \—* B is a surjection from 
£t)[XT] onto E. Hence it extends to an i?-linear map : R[X] —» E. 

Next, by factoring B into a product of primes, we can find n € Z such that B\n. 
Thus, Bl = n for some I G E. Write I = p(B) for some polynomial p in .Eo[X]; we 
thus get that n = Bl = Bp{B) = c m B m+1 + • • • + cqB for some q € Eq, and c m ^ 0. 
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Thus c m is a unit, and we get that 

q(B) = c m B m+1 + • • • + c B - n = 0. 

This means that c^q(X) is in the kernel of the surjection R[X] — > E (which is 
the evaluation map at B). Since c~ 1 g(X) is monic, hence 1, B, B 2 , . . . , B m generate 
E over R, and we are done. □ 

Remarks: The result implies that E is not a polynomial ring -F[-X"]. Moreover, 
if E x is finite, then E is a finitely generated abelian group. Also, the proof only 
assumes that the Euclidean algorithm terminates for all (s, B), where B is any fixed 
element not in Eq, of least v- value. 

Similar questions have been answered; we give a couple of examples. Given n G N, 
we define a Euclidean domain to be of type n if, for every s,d/0 (as in the definition 
of a Euclidean domain E), such that d \ s, there exist exactly n distinct pairs (qi, r,) 
satisfying the Euclidean algorithm property s = qid + r, VI < i < n. We then have 
the following two results. 



Theorem 5 ( [Gap . The only Euclidean domain of type 2 is E = Z. 

Theorem 6 ( RhllJoj). The only Euclidean domains of type 1 are E = F or F[X], 
where F is a field and X is transcendental over F. (In particular, E x = F x .) 

Similar results were also shown by Jacobson and Picavet, in jjaj and [Pij respectively. 



Part 3 : Quantitative aspects (written with P. Moree) 

According to |Kvj . many people, including Johann Bernoulli III, C.-F. Gauss, 
A.H. Beiler, S. Yates, and others, have worked on the problem of repeating deci- 
mals. Repeating decimals were quite a popular topic of study in the 19th century 
(cf. Zentralblatt). We also find the following remarks: 

All prime numbers coprime to 10 can be divided into three groups: 

(1) p such that Dp = p — 1 (the full-period primes). 

(2) p such that D p < p — 1 is odd (the odd-period primes). 

(3) p such that D p < p — 1 is even (the non-full-period primes). 

It was found that the proportion of primes in these three groups had a rela- 
tively stable asymptotic ratio of around 9:8:7, by numerical computations up to 
p = 1370471. We will now explain that it is possible to be more precise than this. 

Concerning the primes in (1). The full-period primes have been well stud- 
ied, starting with C.-F. Gauss in his 1801 masterpiece Disquisitiones Arithmeticae. 
Though no written source for this seems to be available, folklore has it that Gauss 
conjectured that there are infinitely many full-period primes. 
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In September 1927, Emil Artin made a conjecture which implies that the pro- 
portion of primes, modulo which 10 is a primitive root (i.e. primes of type (1) 
in the classification above), should equal a number we now call Artin 's constant 
A = 0.3739558136192 .... More precisely, we have 



J- J- V ntrt - 



Pip ~ 1) 
f 

where the product is over all primes. Clearly, this is close to 9/24 = 3/8 = 0.375, 
the ratio mentioned above. Assuming the Generalised Riemann Hypothesis (GRH) 
it was proved by Hooley |Hoj that the full-period primes have proportion A. 

Artin's original conjecture gave a prediction for the density of primes p such that 
a prescribed integer g is a primitive root modulo p. Following numerical calculations 
by Derrick H. Lehmer and Emma Lehmer in 1957, Artin corrected his conjecture 
for certain g (see |St|). and his corrected conjecture, also attributed to Heilbronn, 
was proved by Hooley in the paper cited above, on assuming GRH. Further general- 
isations of Artin's primitive root conjecture are discussed, e.g. in the survey paper 
[Mo3] . 



In He , Heath-Brown, improving on earlier work by Gupta and Ram Murty GR , 
proved a result which implies, unconditionally, that there exist at most two primes 
(p > 0) - and three squarefree integers k > 1 - for which there are only finitely many 
full-period primes (i.e. for which the qualitative version of Artin's primitive root 
conjecture fails). 

Concerning the primes in (2). Since for an odd prime p, the number p — 1 is 
even, the primes p in (2) can be alternatively described as the primes p such that 
Dp is odd. These primes also have been the subject of study (by mathematicians 
including Sierpinski, Hasse and Odoni). Without assuming any hypothesis, it can 
be shown that the proportion of primes p such that the order of g modulo p with g 
any prescribed integer is odd, exists and is a computable rational number. It turns 
out (cf. |Ba| Theorem 3.1.3] or |Mol| Corollary 1]), that if the base B is not of 
the form ±u 2 or ±2u 2 for any ueN, then the proportion of primes p (among all 
primes) with odd period, is 1/3. Thus the proportion of the primes in (2) equals 1/3. 

Concerning the primes in (3). The proportion of primes such that D p is even 
equals 1 — 1/3 = 2/3. Thus, on assuming GRH, the proportion of primes (3) equals 
2/3 — A. As shown in Lemma |H1 above, B has even period mod p with p \ 2B iff p 
divides B + 1 for some I > 1. Thus the even period condition modulo p is related 
to the divisibility of certain sequences by p. 

In brief, we have arrived at the following result: 

Proposition 4. Assume GRH. Then the natural densities of the sets of primes (1), 
(2), respectively (3), are given in the table below. 
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(1) 


(2) 


(3) 


5 


A 


1/3 


2/3- A 


&d 


0.37395 . . . 


0.33333 . . . 


0.29271... 


Kvant 


9/24 


8/24 


7/24 


/~>-> 


0.37500 . . . 


0.33333 . . . 


0.29166... 



We next look at the residual index c p = [(Z/pZ) x : {B}] (mentioned above), 
which was also equal to the number of distinct chains of recurring decimals (in base 
B). Recall that we had b p c p = p — 1, where b p is the order of B modulo p. 

Note that given a fixed m£N, the number of primes p such that b p = m is finite, 
since there are only finitely many primes less than B m . We now ask for the density 
of the set of primes for which the residual index is fixed, namely {p prime : c p = m}. 
(Note that if m = 1, then this was Artin's conjecture.) 

This too has been answered: in |Mulj we find that assuming the GRH, the den- 
sity of this set Ng equals Cg for certain B, which has been explicitly mentioned 
therein. This density roughly decreases in the order of m~ 2 . Similar results were 
obtained for arbitrary B by Wagstaff |Waj . 

Finally, we ask how often the order of B modulo p differs from that modulo p 2 . 
Related to this is the question: how often does p 2 divide B p_1 — 1 (or not) ? This 
is related to the Wieferich criterion, which first arose in the study of Fermat's last 
theorem. 

For x € N, let us denote the primes less than x, in each of the two sets above, 
by Si(x) and 6*2(2;), respectively. Note that every prime occurs either in Si or in 
5*2. From Mu2j we see, roughly speaking, that the size of 6*2(2;) can be normally 
approximated (if B > 1) by a constant times lnlnx. 

Moreover, if we assume the a&c-conjecture to be true, then we infer (cf. fSIj) that, 
as x tends to infinity, the cardinality of {p < x : p 2 \ (i? p_1 — 1)} exceeds cb In x, 
where cb depends on B but not on x. 

Acknowledgments : Thanks are due to Travis Schedler for finding the counterex- 
ample p = 487 in Section |S] above, to Richard Cudney for first telling me about 
Artin's conjecture, to Professor Dan Shapiro, for telling me about the various ref- 
erences mentioned in Section |H1 above, and to the referee for his suggestions. I also 
especially thank Professor Pieter Moree for his remarks and suggestions, that helped 
bring this paper into its present form. 
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